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We discuss the one-body potentials for protons and neutrons obtained 
from Dirac-Brueckner-Hartree-Fock calculations of neutron-rich matter, in 
particular their dependence upon the degree of proton/neutron asymmetry. 
The closely related symmetry potential is compared with empirical informa- 
<^ ■ tion from the isovector component of the nuclear optical potential. 
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I. INTRODUCTION 



fN| ■ Lately we have been concerned with probing the behavior of the isospin-asymmetric 

equation of state (EOS) [1]. With our recent work on neutron radii and neutron skins [2], 
if} ', we have explored applications of the EOS at densities typical for normal nuclei. 

It is also important to look into systems that are likely to constrain the behavior of 
the EOS at higher densities, where the largest model dependence is observed. Supernova 
explosions and neutron star formation/stability are phenomena where the nuclear EOS plays 
a crucial role. The symmetry energy determines the proton fraction in neutron stars in 
(3 equilibrium, and, in turn, the cooling rate and neutrino emission. Models of prompt 
supernova explosion and systematic analyses of neutron star masses provide often conflicting 
information on the "softness" of the EOS and its incompressibility at equilibrium. 

On the other hand, collisions of neutron-rich nuclei, which are the purpose of the Rare 
Isotope Accelerator (RIA), provide a unique opportunity to obtain terrestrial data suitable 
for constraining the properties of dense and highly asymmetric matter. Such reactions are 
capable of producing extended regions of space/time where both the total nucleon den- 
sity and the neutron/proton asymmetry are large. Isospin-dependent Boltzmann-Uehling- 
Uhlenbeck (BUU) transport models [3] include isospin-sensitive collision dynamics through 
the elementary pp, nn, and np cross sections and the mean field. The latter is a crucial 
isospin-dependent mechanism, and is the focal point of this paper. The contribution to the 
mean field from the neutron/proton asymmetry can be measured through isospin-sensitive 
observables [4]. In summary, this is a timely and exciting topic, which is stimulating new 
effort, on both the experimental and the theoretical sides. 

At this time, it is fair to say that the model dependence of the isospin asymmetric EOS 
is rather large. In fact, even the qualitative behavior of some predictions is controversial, 
as is the case, for instance, with the density dependence of the symmetry energy, upon 
which isospin diffusion in heavy-ion collisions is found to depend sensitively [5]. Thus any 
additional constraint is desirable and should be fully explored. As discussed in Ref. [6], 
nucleon-nucleus optical potential information can be exploited to constrain the strength and 
the energy dependence of the single- neutron/proton potentials in asymmetric nucler matter. 
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The basic idea is that, even though infinite nuclear matter is an idealized system, the single- 
nucleon potentials should bear a clear signature of the optical potential in the interior of the 
nucleus. 

In this paper we will discuss the predictions for the single-neutron/proton potentials and 
the closely related symmetry potential as obtained from our Dirac-Brueckner-Hartree-Fock 
calculations of asymmetric matter. We will compare with other predictions from the liter- 
ature as well as empirical optical potential information. We will point out the large model 
dependence of predictions for those observables that depend sensitively on the difference 
between neutron and proton properties in asymmetric matter. Additional experimental 
constraints are therefore important. Moreover, microscopic, parameter- free approaches are 
the best way to gain deeper insight into the isospin-dependent properties of nuclear matter. 

II. THE SINGLE-NUCLEON POTENTIALS 
A. Momentum dependence 

Unless otherwise specified, we use the Bonn-B potential [7] and the relativistic Brueckner- 
Hartree-Fock (DBHF) model outlined in Ref. [1]. 

We begin by examining the momentum dependence of U n / P , the single neutron/proton 
potential in neutron-rich matter. In Fig. 1, we show U n / p as a function of the momentum 
and for different values of the asymmetry parameter, a = (p n — p p )/(p n + P P ), with p n and 
p p the neutron and proton densities. The total nucleon density considered in the figure is 
equal to 0.185 fm~ 3 and corresponds to a Fermi momentum of 1.4 fm _1 , which is very close 
to our predicted saturation density. 
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FIG. 1. The single-neutron (upper panel) and single-proton (lower panel) potential as a function 
of the nucleon momentum for three different values of the asymmetry parameter. The average Fermi 
momentum is 1.4 fm _1 . 
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For increasing values of a, the proton potential becomes increasingly attractive while 
the opposite tendency is observed in U n . This reflects the fact that the proton- neutron 
interaction, the one predominantly felt by the single proton as the proton density is depleted, 
is more attractive than the one between identical nucleons. Also, as it appears reasonable, 
the dependence on a becomes weaker at larger momenta. 
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FIG. 2. Comparison between DBHF and BHF predictions of the single-neutron (upper panel) 
and single-proton (lower panel) potential. The value of the asymmetry parameter is fixed to 0.4 
and the average Fermi momentum is 1.4 fm _1 . 

In Figure 2 we show the DBHF results in comparison with those from (non-relativistic) 
conventional Brueckner-Hartree-Fock (BHF) calculations. We make the comparison to show 
the considerable difference between the two sets of results as well as to check that our 
BHF predictions are in qualitative agreement with other studies based on the conventional 
Brueckner G-matrix approach. An older work based on that approach can be found, for 
instance, in Ref. [8], where separable representations of the nucleon-nucleon interaction are 
adopted. More recent calculations have been reported in Ref. [9], where the CD-Bonn 
potential [10] is used in conjunction with the BHF approximation. 

The role of the momentum dependence of the symmetry potential in heavy-ion collisions 
was recently examined [11] and found to be important. Symmetry potentials with and 
without momentum dependence and yielding similar predictions for the symmetry energy 
can lead to significantly different predictions of collision observables [11]. 

B. Asymmetry dependence and the symmetry potential 

Regarding U n / P as functions of the asymmetry parameter a, one can easily verify that 
the following approximate relation applies 
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U n / P (k, k F , a) fa U n/p (k, k F , a = 0) ± U sym (k, k F )a (1) 

with the ± referring to neutron/proton, respectively. Figure 3 displays the left-hand side of 
Eq. (1) for fixed density and nucleon momentum and clearly reveals the linear behaviour of 
U n / P as a function of a. 




FIG. 3. The single-neutron (upper panel) and single-proton (lower panel) potential as a function 
of the asymmetry parameter for fixed average density (k F = 1.4 fm -1 ) and nucleon momentum (k— 
k F ). 



Although the main focus of Fig. 3 is the a dependence, predictions are displayed for 
the Bonn A, B, and C potentials [7]. These three models differ mainly in the strength of 
the tensor force, which is mostly carried by partial waves with isospin equal to and thus 
should fade away in the single-neutron potential as the neutron fraction increases. Reduced 
differences among the three models are in fact observed in U n at the larger values of a. 

Already several decades ago, it was pointed out that the real part of the nuclear optical 
potential depends on the asymmetry parameter as in Eq. (1) [20]. Thus, the quantity 



U n + Up 



= U , (2) 



which is obviously the single-nucleon potential in absence of asymmetry, should be a reason- 
able approximation to the isoscalar part of the optical potential. The momentum dependence 
of Uq (which is shown in Fig. 1 as the a=0 curve), is important for extracting information 
about the symmetric matter EOS and is reasonably agreed upon [5,12-19]. 
On the other hand, 

Un ~ U P - TJ (*\ 
2 a ~ u sym {O) 
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should be comparable with the Lane potential [20], or the isovector part of the nuclear op- 
tical potential [20]. (Notice that in the two equations above the dependence upon density, 
momentum, and asymmetry has been suppressed for simplicity.) We have calculated U sym 
close to nuclear matter density and as a function of the momentum, or rather the corre- 
sponding kinetic energy. The predictions obtained with Bonn A, B, and C are shown in 
Fig. 4. They are compared with the phenomenological expression [20] 

U L ane = CL - bT (4) 

where T is the kinetic energy, a w 22 - 3AMeV, b « 0.1 - 0.2MeV. 




FIG. 4. The symmetry potential as a function of the nucleon kinetic energy at nuclear matter 
density. The predictions obtained with Bonn A, B, and C are compared with empirical information 
from nuclear optical potential data (shaded area). See text for details. 



The strength of the predicted symmetry potential decreases with energy, a behavior 
which is consistent with the empirical information. The same comparison is done in Ref. [6] 
starting from a phenomenological formalism for the single-nucleon potential [21,22]. There, 
it is shown that it is possible to choose two sets of parameters which lead to similar values 
of the symmetry energy but exactly opposite tendencies in the energy dependence of the 
symmetry potential as well as opposite sign of the proton-neutron mass splitting. As a 
consequence of that, these two sets of parameters lead to very different predictions for 
observables in heavy-ion collisions induced by neutron-rich nuclei [22]. 

Our effective masses for proton and neutron are shown in Fig. 5 as a function of a 
and at saturation density. The predicted effective mass of the neutron being larger than 
the proton's is a trend shared with microscopic non-relativistic calculations [8]. In the non- 
relativistic case, one can show from very elementary arguments based on the curvature of the 
single-particle potential that a more attractive potential, as the one of the proton, leads to a 
smaller effective mass. In our DBHF effective-mass approximation, we assume momentum- 
independent nucleon self-energies, Us and Uy, with a vanishing spacial component of the 
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vector part. In such limit, following similar calculations of symmetric matter [23], the one- 
body potential is written as [1] 

771* 

Ui(p) = -j£Us,i + U v ,i (5) 



where E* = y (m*) 2 + p 2 , m* = rrii+Us,i, and i = n or p for neutrons or protons, respectively. 
Defining for convenience U 0ti = Ug,i + Uy,i, the expression above becomes a two-parameter 
formula which requires the fitting of two constants, just like in the non-relativistic case. 
Now, since the single-proton potential is more attractive (see Fig. 1), and both the neutron 
and proton potentials tend to the same limit at high momenta, it is easy to see from Eq. (5), 
or rather its derivative, that the proton effective mass obtained in this way must be smaller 
than the neutron's. 




FIG. 5. The proton and neutron effective mass as a function of the asymmetry parameter and 
for fixed average density (&f= 1-4 fm _1 ). 



One can encounter in the literature the statement that relativistic mean field (RMF) 
models predict a switch in sign of the neutron-proton effective mass splitting whenever the 
scalar isovector 5 meson is included [22]. We suggest that the reason for this mechanism is 
in the first-order nature of RMF calculations. The first-order contribution of the S meson 
is known to be repulsive in 3 Si (the most important partial wave for generating nuclear 
binding). Therefore, in a context where the proton density increases (thus making the np 
interaction the dominant contribution to the single-proton potential), a first-order calcula- 
tion would in fact generate additional repulsion in the single-proton potential. Hence the 
larger proton effective mass. On the other hand, in more microscopic approaches, such as 
the one we use, the nucleon potential is iterated to all orders via the in-medium scattering 
equation. In particular, most of the attraction in the np interaction is generated in second 
order, through the two-pion exchange. (The 5 meson is included, of course, but mainly 
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for the sake of accurate reproduction of NN phase shifts.) As we already observed in Sec- 
tion IIA, increasing the np contribution to the single-proton potential results in increased 
attraction. 




FIG. 6. Contribution from the asymmetry to the average potential energy per neutron (upper 
panel) and proton (lower panel). Average density as in the previous figures. 



Some comments are in place concerning other DBHF calculations of asymmetric matter 
where the Dirac neutron effective mass is reported to be larger than the proton's [25,26]. In 
Ref. [25], the authors argue that the Dirac mass, m* = to + Us, should not be compared 
with the effective mass which can be extracted, for instance, from analyses based on non- 
relativistic optical models. Instead, an effective mass based on the energy dependence of the 
Shroedinger equivalent potential should be calculated, in which case one finds that m* > to*. 
The arguments found in Ref. [25] are well known and were already advanced in Ref. [27], 
where as many as six different definitions for the effective mass are introduced. Again, one 
must keep in mind the point we already made above. If the nucleon self-energy is written in 
terms of a scalar potential and only the time-like component of a vector potential, and both 
are taken to be momentum independent, it is then easy to see that the expansion of the 
single-particle energy is consistent to leading order with the non-relativistic single-particle 
energy. Thus it is reasonable, and in fact to be expected, that our Dirac masses would be 
qualitatively consistent with those from non-relativistic predictions, such as BHF calcula- 
tions. In summary, the effective mass is just part of a convenient parametrization of the 
single-particle potential. Clearly, how many terms are retained in the nucleon self-energy 
and how their momentum dependence is handled will impact the parametrization. Ulti- 
mately, physical observables depending on the neutron/proton mean field must be correctly 
described, irrespective of the chosen parametrization. 

One conclusion that can be made from all of the above is that the parameters of the single- 
nucleon potential in asymmetric matter appear to be weakly correlated to observables such 
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as the energy per particle or the symmetry energy, where proton and neutron contributions 
are averaged together. Constraints from "differential" or relative observables, namely those 
specifically sensitive to the difference between proton and neutron properties, are thus very 
much needed. 

Before closing, we also show for completeness the average potential energy per neu- 
tron/proton, where the momentum dependence has been integrated out. This is the pro- 
ton/neutron potential energy contribution to the total energy per particle which then ap- 
pears in the EOS. Actually, what we show in Fig. 5 are the average potential energies from 
which the part coming from the symmetric EOS has been subtracted out, that is, just the 
contribution from the asymmetry to the interaction potential energy, 

< AU n/p > (p, a) =< U n/p > (p, a)- < U(p, a = 0) > . (6) 

Clearly, the contribution from the asymmetry, in both the momentum-dependent and the 
momentum-averaged potentials, turns out to be large and positive for neutrons, large and 
negative for protons. This component of the mean field will then be effective in separating 
the collision dynamics for neutrons and protons by making more neutrons unbound than 
protons (or, by making the neutrons more energetic, if already unbound). This effect can 
be discerned through observables such as the neutron/proton differential flow in heavy-ion 
collisions [4]. 

III. CONCLUSIONS 

We have focussed on some of the properties of neutrons and protons in neutron-rich 
matter. This is a topic of contemporary interest. Its relevance extends from the dynamics 
of colliding nuclei to nuclear astrophysics. 

Different models may be in fair agreement with respect to averaged properties of the EOS, 
and yet produce very different predictions of properties such as the symmetry potential and 
the closely related single-nucleon potentials and effective masses. Clearly, more stringent 
constraints are needed for the isospin-dependent properties of the EOS. 

Very good transport model calculations are available from the literature [3,4]. However, 
considerable amount of phenomenology is often involved in the input of these models (for 
instance, the mean field is based on some phenomenological interaction [24,11] and/or the 
elementary cross sections are obtained from empirical data). We calculate all of the above 
ingredients microscopically and internally consistent with respect to the two-body force. 
We are presently studying the dependence on density and asymmetry of the in-medium 
isospin-dependent nucleon-nucleon cross sections with the purpose of obtaining a convenient 
parametrization as a function of energy, density, and asymmetry. Our microscopic informa- 
tion (both elementary cross sections and mean field) can be a valuable input for transport 
model calculations of heavy-ion dynamical observables. This combined effort will comple- 
ment new data to be taken at RIA and eventually shed light on the less known aspects of 
the nuclear equation of state. 



ACKNOWLEDGMENTS 



8 



The authors acknowledge financial support from the U.S. Department of Energy under grant 
number DE-FG02-03ER41270. 



D. Alonso and F. Sammarruca, Phys. Rev. C 67, 054301 (2003). 

D. Alonso and F. Sammarruca, Phys. Rev. C 68, 054305 (2003). 

See, for instance, B.A. Li, Phys. Rev. Lett. 85, 4221 (2000), and references therein. 

B.A. Li, CM. Ko, and Z. Ren, Phys. Rev. Lett. 78, 1644 (1997); Phys. Rev. Lett. 88, 192701 

(2002). 

L.-W. Chen, CM. Ko, and B.-A. Li, nucl-th/0407032. 

B. A. Li, Phys. Rev. C 69, 064602 (2004). 

R. Machleidt, Adv. Nucl. Phys. 19, 189 (1989). 

I. Bombaci and U. Lombardo, Phys. Rev. C 44, 1892 (1991). 

Kh.S.A. Hassaneen and H. Miither, nucl-th/0408035. 

R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

Bao-An Li, Champak B. Das, Subal Das Gupta, and Charles Gale, Phys. Rev. C 69, 011603 
(2004). 

C. Gale, G. Bertsch, and S. Das Gupta, Phys. Rev. C 35, 1666 (1987). 

CM. Welke, M. Prakash, T.T.S. Kuo, S. Das Gupta, and C. Gale, Phys. Rev. C 38, 2101 
(1988). 

C. Gale, CM. Welke, M. Prakash, S.J. Lee, and S. Das Gupta, Phys. Rev. C 41, 1545 (1990). 
Q. Pan and P. Danielewicz, Phys. Rev. Lett. 70, 2062 (1993). 

J. Zhang, S. Das Gupta, and C. Gale, Phys. Rev. C 50, 1617 (1994). 

V. Greco, A. Guarnera, M. Colonna, and M. Di Toro, Phys. Rev. C 59, 810 (1999). 

P. Danielewicz, Nucl. Phys. A673, 375 (2000). 

D. Persram and C. Gale, Phys. Rev. C 65, 064611 (2002). 

A. M. Lane, Nucl. Phys, 35, 676 (1962). 

I. Bombaci, Chap. 2 in Isospin Physics in Heavy-Ion Collisions at Intermediate Energies, Eds. 

B. A. Li and W. Udo Schroder (Nova Science Publishers, Inc., New York, 2001). 
J. Rizzo et al, Nucl. Phys, A732, 202 (2004). 

R. Brockmann and R. Machleidt, Phys. Lett. B 149, 283 (1984); Phys. Rev. C 42, 1965 (1990). 

C. B. Das, S. Das Gupta, C. Gale, and Bao-An Li, Phys. Rev. C 67, 034611 (2003). 
Z-Y. Ma, J. Rong, B.-Q. Chen, Z.-Y. Zhu, H.-Q. Song, Phys. Lett. B, in press. 

E. N.E. van Dalen, C. Fuchs, and A. Faessler, nucl-th/0407070. 
M. Jaminon and C. Mahaux, Phys. Rev. C 40, 354 (1989). 



9 



